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Abstract — In this article we obtain estimates on tlie approx- 
imate eigenstructure of cliannels with a spreading function 
supported only on a set of finite measure \U\. Because in typical 
application like wireless communication the spreading function 
is a random process corresponding to a random Hilbert-Schmidt 
channel operator Ti. we measure this approximation in terms of 
the ratio of the p-norm of the deviation from variants of the 
Weyl symbol calculus to the a-norm of the spreading function 
itself. This generalizes recent results obtained for the case p = 2 
and a = 1. We provide a general approach to this topic and 
consider then operators with \U\ < oo in more detail. We show 
the relation to pulse shaping and weighted norms of ambiguity 
functions. Finally we derive several necessary conditions on |L'^|, 
such that the approximation error is below certain levels. 

I. Introduction 

Optimal signaling through linear time-varying (LTV) chan- 
nels is a challenging task for future communication systems. 
For a particular realization of the time-varying channel opera- 
tor the transmitter and receiver design which avoids crosstalk 
between different time-frequency slots is related to "eigen- 
signaling". Eigen-signaling simplifies much the information 
theoretic treatment of communication in dispersive channels. 
However, it is well-known that for a whole class of channels 
such a joint separation of the subchannels can not be achieved. 
A typical scenario, present for example in wireless com- 
munication, is signaling through a random doubly-dispersive 
channel Ti.: 

r{t) = {ns){t)+n{t) 

From signal processing point of view the preferred design of 
the transmit signal s{t) needs knowledge on the true eigen- 
structure of Ti.. This would in principle allow interference- 
free transmission and simple recovering algorithms of the 
information from received signal r{t) degraded by the noise 
process n{t). However, for Ti. being random, random eigen- 
structure has to be expected in general and a joint design of the 
transmitter and the receiver for an ensemble of channels has to 
be performed. Nevertheless with such an approach interference 
can not be avoided and remains in the communication chain. 
For such interference scenarios it is important to have bounds 
on the distortion of a particular selected signaling scheme. 

First results in this field can be found already in the literature 
on pseudo-differential operators [1], [2]. More recent results 
with direct application to time-varying channels were obtained 
by Kozek [3] and Matz [4] which resemble the notion of 
underspread channels. They investigated the approximate sym- 
bol calculus of pseudo-differential operators in this context 



and derived bounds for the ^2-norm of the distortion which 
follow from the approximate product rule in terms of Weyl 
symbols. Controlling this approximation intimately scales with 
the "size" of the spreading of the contributing channel oper- 
ators. For operators with compactly supported spreading this 
is \U\ - the size of the spreading support U. Interestingly this 
approximation behavior breaks down in their framework at a 
certain critical size. Channels below this critical size are called 
in their terminology underspread, otherwise overspread. 

Underspreadness of time-varying channels occurs also in 
the context of channel measurement [6]. See also the recent 
article [7] for a rigorous treatment of channel identification 
based on Gabor (Weyl-Heisenberg) frame theory. The authors 
connect the critical time-frequency sampling density imma- 
nent in this theory to the stability of the channel measurement. 
A relation between these different notions of underspreadness 
has to be expected but will be out of the scope of this paper 

This article considers the problem of approximate eigen- 
structure from a different angle, namely investigating the £p- 
norm Ep of the error Ti.s — \r for well-known choices of A. 
This direct formulation allows for improvements to the exist- 
ing bounds, generalizations to arbitrary regions of spreading 
and different distortion measures. Furthermore the approach 
will show the connection to well-known fidelity criteria related 
to pulse design [3], [8], [9]. Using the techniques recently 
presented in [9] we finally extract necessary conditions for \U\ 
that the error does not exceed certain levels. Furthermore we 
discover some interesting relations related to underspreadness 
in form of [3], [4]. 

The paper is organized as follows. After settling the basic 
definitions in the first section of the paper the second section 
reviews the Weyl correspondence and the spreading represen- 
tation of Hilbert-Schmidt operators. In the next sections we 
will then consider the problem of controlling Ep and give 
necessary conditions on \U\. In the last part we will verify 
our framework with some numerical tests. 



A. Some Definitions 

But before starting, the following definitions are needed. 

def 



For 1 < p < oo and / : R — > C the functional || 
(/ \fi't)\'^dtY^^ is then usual notion of the p-norm (dt is the 

def 

Lebesgue measure on R). Furthermore for p = cxa is ||/||oo = 
ess sup \f{t)\. If ll/llp is finite / is said to be in £p(R). We 
will frequently make use of the relation = ||/||ac- The 



function xu will always denote the characteristic function onto 
the set U CR'^. 

II. Displacements Operators and Ambiguity 
Functions 

Time-frequency representations are an important tool in sig- 
nal analysis and physics. Among them are Woodward's cross 
ambiguity function and the Wigner distribution. Ambiguity 
functions can be understood as inner products representations 
of displacement (or shift) operators, defined by its action on 
function / : M ^ C as; 

iSJ)ix)■.= e^^-^-^f{x~fI^) (1) 

The functions / are time-signals, i.e. ^ ~ {111,^2) G 
where /^i is time shift and 112 is a frequency shift. However the 
following can be straightforward extended to multiple dimen- 
sions. The operator acts isometrically on all £p(K), hence 
is unitary on the Hilbert space £2(K)- The arbitrariness due to 
the non-commutativity of shifts in the previous definition can 
be covered in using generalized displacements of the form: 

S^{a) = S(o,^2(1+q))S'(^j,o)S'(o,p2(1_q)) (2) 

We will call a as polarization. All these operators establish (up 
to unitary equivalence) unitary representations of the Weyl- 
Heisenberg group on £2(1^) (see for example [2]). In physics 
it is common to choose the most symmetric case a = and 
the operators are usually called Weyl operators. The definition 
O appears for a = 1/2 and we call 5'^ = 5'^ (1/2) as time- 
frequency shift operator. If we define the symplectic form as 
r]{^,v) := ^iiV2 — ^2Vi, we have the following well-known 
Weyl commutation relation: 

S^{a)S,{(3) - e-*2.„(M,-)5^(^)5^(^) (3) 
In this way the cross ambiguity function can be defined as: 

III. Weyl Correspondence and the Spreading 
Representation 

In the previous section it was motivated that Ag"-* yields 
a local time-frequency description of functions. Now the 
same can be repeated for Hilbert-Schmidt operators. Let us 
introduce them as the 2th Schatten class: If we define for 
linear mappings A £ L(£2(1R)) from Hilbert space >C2(K) 
into itself \A\ := {A*A)'^/^, then for 1 < p < cx) the 
functional \\A\\p := Tr(|A|^')"^^^ is called the pih Schatten 
norm. The set Tp := {A e L{C2{M.)), \\A\\p < 00} is called 
the pth Schatten class where 7^ is set to be the compact 
operators. Then Tp for 1 < p < 00 are Banach spaces and 
%_ C Tp C Too (see for example [10] or [11]). The sets Ti 
and T2 are called trace class and Hilbert-Schmidt operators. 
Hilbert-Schmidt operators form itself a Hilbert space with 
inner product {A,B)t2 = TrA*B. 

In particular for any Ti. <E Ti there holds by properties of the 
trace Tr(X7i) < ||?<||i||X||, where ||-|| denotes the operator 
norm. Hence for X = Sn{a) given by (|2]i one can define with 



analogy to ordinary Fourier transform [12], [13] a mapping 

F(") : Ti ^ C2{R'^) via 

(F(")7i)(/i) Tr{S;{a)n) = {S,{a),n)r, (4) 

Note that (F("'7i)(0) = TrH and |(F("''H)(/^)| < \\n\\i. 
The function F'"^7i G £2(K^) is sometimes called the "non- 
commutative" Fourier transform [11], inverse Weyl transform 
[14] or a-generalized spreading function of Ti. [3]. 

Lemma 1 (Spreading Representation) Let Ti e 7^. Then 
there holds 

1 (F(")7i)(/i)5^(a)d// = / {S^ia),n)T,S^{a)dfM 

(5) 

where the integral is meant in the weak sense. 

The extension to 72 is due to continuity of F^") : 7i 
£2(R^) and density of 7i in 72. A complete proof of this 
lemma can be found for example in [11]. For Ti, X G 72 the 
following Parseval-like identity 

{x,n)T, = {F'-"^x,F'-"^n) (6) 

holds. If we define the symplectic Fourier transform of a 
function JF^ : ^ C as: 

{TsF){fi) = [ e-'2-''(-.A')F(i.)di. (7) 

then JTgF'"' establishes a correspondence between the ordi- 
nary function L^^'' = JRjF^^^X and an operator X (Weyl 
quantization [14]). The function i^"* is called (generalized) 
Weyl symbol of X. The original Weyl symbol is L^x^ . The 
cases a = 1/2 and a = —1/2 are also known as Kohn- 
Nirenberg symbol (or Zadeh's time-varying transfer function) 
and Bello's frequency-dependent modulation function [5]. 
Using Parseval identity for eq. (|6]l extends now to 

{X,Y)t, = (F(")x,F(")r) = {l'^x^I^y^) (8) 

and consequentially ||X||2 = ||F(")X||2 = \\lP\\2. 

iv. elgenstructure of operators with compactly 
Supported Spreading 

A. The Approximate Eigenstructure 

It is of general importance how much the Weyl sym- 
bol or a smoothed version of it approaches the eigenvalue 
characteristics of a given Hilbert-Schmidt operator. Since 
G 72 is compact, it has a Schmidt representatiorQ Ti = 
Hn=i ■'^k{xk,-)yk with the singular values {sfc} of Ti and 
orthonormal bases {xk] and {uk}- However, the latter depends 
explicitely on Ti and can be very unstructured. We are 
interested in a choice which is "more independent" of Ti, 
which give rise to the following definition of what we will 
call "approximative" eigenstructure: 

'For T-L given in matrix representation also known as "Singular Value 
Decomposition" 



Definition 2 (Approximate Eigenstructure) Let \ E C If 

for 5,7 : R — > C with \\g\\p = ||7||p = 1 holds 

||7i7-A.9||p<e(A,7,<?) (9) 

we call A an "£p-approximate eigenvalue " of H with bound 
e(A,7,g). 

Obviously we have e{sk, Xk, Vk) = for each p. The question 
is, how much particular choices for A(^) and functions 
and Sf^g can approach the eigenstructure of a Hilbert Schmidt 
operator 7t. Because in general these functions differ from the 
Schmidt representation of Ti, they will give an error which we 
have to control. Summarizing: 



C2: "B = 1", such that X{fi) ^ This case is related 

to the symbol calculus and needed for comparisons with 
the previous results given so far 
Then the following theorem parallels Theorem [3] and its 

consequence ( fT2] i. 

Theorem 4 For 1 < p < 00 and 1 < a < 00 holds: 



EP 



^ < pP^mi + \B\^ 2/?.{A(")B})xr/|lvp (14) 



where l/a + l/b ~ 1. The minimum over B is achieved for 
B = A ' 



(") 



Proof: The proof follows from the middle term of ( [23] l in 
Lemma |6] given in the next section if one set W ~ and 



Ep-.-WHS^j \{fi)Spg\\p<e{X{^i),S^j,Sf,g) (10) ^ The constant poo will also be explained later on 



At this point we introduce furthermore the following abbrevi- 
ation: := F^"'?^ will always be the spreading function 
of 7i. With ?7 C we will denote its support and with 
\U\ ~ WxuWi its size (measure). 

B. Previous Results 

In [3, Theorem 5.6] W. Kozek has been considered the case 
A = = Ts'S'^ and g = 7. He obtained the following 
result: 



Theorem 3 (W. Kozek [3]) Let U — [—tq^tq] x [-z/Oji^o]. If 
\U\ = 4x01^0 < 1 then 



El<2 sin( 



^|C/|, 



where = ||(A^°^^ - l)xc/||oo- 



■'■H 111 



will 



+ 2||S, 



(0)||2 



■H 111 



(11) 



Using < llS^^llf eq. O can be written as 



.(0)||2 



El 



(0)„2^2sin(^) + 3||(A(°)-l)xd 



(12) 



■'■H 111 



G. Matz generalized the result of Theorem |3] in [4, Theorem 
2.22] to a formulation in terms of weighted 1-moments of 
spreading functions which includes now different polarizations 
a and is not restricted to the special choice of U. For U = 
[—To, To] X [— i^oi t'o] ™d a = the bounds agree with (fT2l i. 

C. New Related Results 

Instead of directly considering the case of using the Weyl 
symbol for the approximation of the eigenstructure, we use a 
"smoothed" version: 



(13) 



and consider two cases 
CI: "B 



Ag"'", such that A = i^'' * JF^A^"' where * 
denotes convolution. This corresponds to the well-known 



smoothing with the cross Wigner function JF^Ag"'. 

^Note that can be replaced with S^(/3) without change of Ep 



It follows that for C2, p = 2 and a = 1 that 

(0<2||(l-A<,^)xc/||oo (15) 

which improves the previous bounds ( fTTT ) and ( fT2] i. It is 
independent of the polarization a and does not require any 
shape or size constraints on U. Interestingly the offset in 
(fT2] |. which does not depend on (.9,7) and in a first attempt 
seems to be related to the notion of underspreadness, has been 
disappeared now. 

V. Generalization and Proofs 

For the study of random operators we have to classify the 
overall spreading function. Thus we assume that all realiza- 
tions of the spreading function can be written as 

^^^\li)^K{ti).W{ti) (16) 

for a common function K : M+, which model 

some apriori knowledge (for example the square root of the 
scattering function in the WSSUS assumption [5] directly 
or some support knowledge). We will always denote with 
J7 C M2 then the support of K. The function W : R'^ ^ C 
represents the random part. From this considerations it is 
desirable to measure the error Ep with respect to a certain 
a-norm ||VF||a of the random part, thus to look at the ratio 
i?p/||M^||a. We have the following Lemma: 

Lemma 5 Let Pp{v) := \\Sv{a)-f - B{v)g\\pK{v). For 1 < 
p<oo, 1 < a < 00 and 1/a + 1/6 = 1 holds 



Epl\\W\\a<\\Pp\\b 

whenever W G Ca{R^^) and Pp G ^^(R^). 



(17) 



Proof: Firstly - using Weyl's commutation rule and defini- 
tion of A in ([13]) gives us 

Ep^wJ di,i:[^\v)SM)St.l-\{p)S^g\\p 



1 II j duT,'-^\v)e-^^-^^^^'^\SMh ^ B{u)g)\\p 



(18) 



Note that p-norm is with respect to the argument of the 
functions g and Sy{a)^. The last step follows because S^{a) 
acts isometrically on all £p(M). If we define 

f{x,u) := e-'2-'(-^)l](?'(i^)[(5,(a)7)(a;) -i?(i^).g(x-)] 

(19) 

eq. (fTSl l reads for 1 < p < cx) by Minkowski (triangle) 
inequality 

Ep^wj dvf{-,v)\\j, < II / d,,\fi;:.)\\\, < ( dv\\f{;v)\\j, 

^^^^ 

With = W ■ K and Holder's inequality follows the claim 
of this lemma. ■ 

Let us fix for the moment II7II2 = II5II2 = 1 (||Ag"''||oo < 1) 
and a constant C such that max(||7||p, \\g\\q) < C. Then it is 
obvious that for CI and C2 follows Pp{i') < 2GK{v) and we 
have: 

i^p/||VK|U < 2C||/i||6 (21) 

In the next Lemma we will show that \\pp\\b can be related to 
weighted norms of ambiguity function, which we have studied 
already in [9]. A central role will play here the function: 



p{v) sup |(S'^(a)7)(a;) - B{v)g{x)\xu{v) 



(22) 



which characterize the relative smoothness of g and 7 with 
respect to shifts ly <E U. Let us furthermore denote its 
supremum with p^o ■= \\p\\oc- Due to limited space we have 
to postpone a detailed discussion of p and poo (which will 
be important for p 2 and then one has also to consider 
Poo = Poo{U)) to a separate journal paper in preparation. 

For simplicity we now make w.l.o.g. the assumption ||g||2 = 
II7II2 = 1 and define the non-negative function R := 1 + 



. (a)|2 



Lemma 6 With the assumptions of Lemma |5] holds: 

p-2 

llPpllb < Poi I 



(23) 



with equality for p — 2. The minimum over B of the rhs is 
achieved for C2. 



i/p 



Proof: We have for p > 1: 

PpH < l^j {\{SM)l){^)-B{v)g{x)\^dxj A» 

= {-p{v)K{v)Y-^ p^{vf'P 

(24) 

with equality for p = 2 such that 

\\Pph < p'S'\\pIK''-X/1 pfWRKX'l (25) 

From the definition of R it is obvious that the minimum of 
the bound in (|24l) is taken at B{U) = aI"\u) which is 
provided by CL Because equality for p = 2 in (l24l i this is 
also the optimizer for ||p2||6 for any b. ■ 



A. Relation to Weighted Ambiguity Norms 

Now we will discuss the connection to weighted norms 
of ambiguity functions and fidelity criteria related to pulse 
shaping as introduced in [9]. It will give (partially) new 
insights into the terms of underspreadness in this context. 
Assume that Roo := ||i?Xf7||oo < 1 and b > p. Then it can be 
shown that 

i/b 



iPpWb 



RK 



(26) 



If the latter can not be fulfilled, hence for 6 < p or if i?oo > 1, 
it still holds: 



\Pp\\b < P 



Rlinm 



(27) 



It can be verified that for CI the condition R^o < 1 is always 
fulfilled. Thus, in this case ( |26l ) reads: 



\Pp\\b < PoS 



|A-^||i-||(A(")f7^^||i 



i/b 



(28) 



We conclude that with our assumptions a maximization of 
the "2-channel fidelity" \\{A[%''fK^\\i (the case r = 2 in 
[9]) controls E'p/ljVFlla. This is also the term important for 
pulse shaping with respect to scattering function of WSSUS 
channels [3], [8]. 

But for the condition C2 the behavior is different. We arrive 
at the very interesting condition, that R^o < 1 is equivalent to 



1 



- < inf Re{A^'^^(A.)} 



(29) 



We will show later on that this condition can not be fulfilled 
on every U . However, if ( |29] | holds we get from (|26] | and 



\\Pp\\b < Po 



2(||/^''||i-||Re{A(^^)}A^''||i 



1/6 



(30) 



which is then a problem of the maximization of the "1-channel 
fidelity" (the case r = 1 in [9]). 



VI. Some Necessary Support Conditions 
In this section we present some necessary condition on |[/| 
for the case K = xu (therefore support conditions on S^') 
which follow from the methods presented in [9]. Due to limited 
space we have to omit the proofs, which will then appear 
separately in a journal version. In particular one can show 
from [9] that for \U\ < emiii(l,2/r) follows 



lA^T) 



-97 I XU\\l 

and with this result follows: 



< \U\e- 



(31) 



Lemma 7 (Necessary Condition for C2) The condition in 
can only be fulfilled if \U\ < 2eln2. 



We believe that this bound is very coarse and it should be 
possible to improve it in using more advanced techniques. 
Further support results follow for b > p. We define from (1281 ) 
the following quantity 

niU) pT {\U\ - ||(A(^)fxc/||i)'^' (32) 
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Fig. 1. Approximation error Ea/IIS^'lb for the case C2. "General 
bound" refers to )2U . "Optimized Bound (Roo = 1)" is )27t and "necessary 
condition" is the results of (34). The monte carlo data is obtained from 
Gaussian signaling and a rectangular' spreading function with independent 
complex normal distributed components Cfe (K = 10) as explained in (35). 

for the case CI. For C2 we have to guarantee that i?oo < 1 
and we define instead from ( [30l ) the quantity: 

r2{U) ■.= pT (2(|C/|-||Re{A(^)}xc/||i))'^'' (33) 
Then the following can be shown: 

Lemma 8 (Necessary Conditions on \U\) Let 1 > Sk > 0, 
\U\ < e and k = {1,2}. If rk {U) < 4 then U has to fulfill: 

P-2 / |!7| \ i/b 

PocT ik\U\{l-e- — )j <4 (34) 

The latter is an implicit inequality for \U\. However, it is 
possible to obtain from this an explicit upper bound for \U\ 
for a given Sk (not shown in the paper). 

VII. NuiviERicAL Verification 

In the following we will evaluate and test the obtained 
bounds for Gaussian signaling, i.e. g and 7 are time-frequency 
symmetric Gaussian functions. We consider a spreading func- 
tion given as: 

EH = ^ CkXQiiy - u{k + o)) (35) 

fcez|- 

where Zk ^ {0 ■ ■ -K-l}, Q = [0, u] x [0, u] and o = (i, i). 
If we fix the support of the spreading function to be \U\, 
then follows u = y/\U\/K. For such a model the a-norm of 
the spreading function is: |lS||a = it^^°||cl|a, where ||c||a is 
simply the ath vector norm of the vector c with coefficients 
Cfc. The error Ep can be simplified much for Gaussian sig- 
naling and finally computed numerically (therefore Gaussian 
signaling was chosen). Fig. [T] shows the results of several 
monte carlo runs, each corresponds to one point in the plot. 
For comparison the various bounds and results are included 



in the plot (more details in the caption). Clearly the most 
important result ("necessary condition") can not serve as a 
bound. However, its interesting that it produce a rough value 
of the approximation error. 

VIII. Conclusions 

In this contribution we established in a more general fashion 
the problem of approximate eigenstructure of LTV channels. 
We extracted several criteria related to signaling in those chan- 
nels and pulse shaping. We hope that our results give some 
more implications to the role of underspreadness for wire- 
less communication. Furthermore the connection to symbol 
calculus of pseudo-differential operators is straightforward, 
such that insights into topics like approximate commutativity 
— or more generally speaking — approximations to spectral 
properties have to be expected. 
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